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Abstract. We study the population of compact X-ray sources in the Large Magellanic Cloud using the archival data of the 
XMM-Newton observatory. The total area of the survey is as 3.8 square degrees with a limiting sensitivity of as 10~ 14 erg/s/cm 2 , 
corresponding to a luminosity of a 3 • 10 33 erg/s at the LMC distance. Out of ~ 460 point sources detected in the 2-8 keV 
energy band, the vast majority are background CXB sources, observed through the LMC. Based on the properties of the optical 
and near-infrared counterparts of the detected sources we identified 9 likely HMXB candidates and 19 sources, whose nature 
is uncertain, thus providing lower and upper limits on the luminosity distribution of HMXBs in the observed part of LMC. 
When considered globally, the bright end of this distribution is consistent within statistical and systematic uncertainties with 
extrapolation of the universal luminosity function of HMXBs. However, there seems to be fewer low luminosity sources, 
log(Lx) S 35.5, than predicted. We consider the impact of the "propeller effect" on the HMXB luminosity distribution and 
show that it can qualitatively explain the observed deficit of low luminosity sources. 

We found significant field-to-field variations in the number of HMXBs across the LMC, which appear to be uncorrelated with 
the star formation rates inferred by the FIR and H a emission. We suggest that these variations are caused by the dependence of 
the HMXB number on the age of the underlying stellar population. Using the existence of large coeval stellar aggregates in the 
LMC, we constrain the number of HMXBs as a function of time t elapsed since the star formation event in the range of r from 
-1-2 Myr to ~ 10 - 12 Myr. 
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1. Introduction. 

First imaging X-ray obs ervations of the La rge Magellanic 
Cloud with Ein s tein .Wang et all Il99ll) and ROSAT 
(Haberl & Pietschl Il999l) observatories revealed a moder- 
ate population of several tens of X-ray sources associated with 
this closest neighbor of the Milky Way. While in the soft X-ray 
band a large fraction of these sources are supernovae remnants, 
at higher energies, above ~ few keV, X-ray binaries provide 
the dominant contribution. For examp le, based on th e spectral 
hardness of bright ROSAT sources, Kahabka (2002) identified 
in the entire LMC ~ 20 - 30 probable X-ray binary candidates 
with luminosity exceeding Lx ^ few x 10 35 erg/s. 

The relative numbers of high and low mass X-ray bina- 
ries (H or LMXBs) are defined by the specif i c star formation 
rate S FR/M„ of the host galaxy (Gri mm et aT_l2003_lGilfanovi 
2004). Owing to rather low mass, M* ~ 2 • 10 9 M© (section 
I3.2I . and moderate star formation rate (SFR) of ~ 0.5 Mq 
yr~' (section l3~3*t . this quantity is rather high for the LMC, 
with a SFR/M, ~ (2 - 3) • 10 -10 yr _1 , exceeding by a fac- 
tor of ~ 5 - 10 that of our Galaxy. Correspondingly, the ratio 



of the expected numbers of bright, log(Lx) ~ 35, X -ray bina- 
ries i n the LMC equals Nlmxb '■ Nhmxb ~ 1 : 7 jGirfanovt 
2004), i.e. is nearly opposite to that observed in the Milky Way 
dGrimm et all l2002ll . The absolute number of HMXBs in the 
LMC is expected to be ~ 1 /4 of the number observed in our 
Galaxy, corresponding to the ratio of star formation rates in 
these two galaxies. 
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As has been shown bv iGrimm et all d2003l) . the X-ray lu- 
minosity function (XLF) of HMXBs obeys, to the first approx- 
imation, a universal power law distribution with a differential 
slope of « 1 .6, whose normalization is proportional to the star 
formation rate of the host galaxy. The validity of this univer- 
sal HMXB XLF has been established in a broad range of the 
star formation rates and regimes and in the luminosity range 
\og(Lx) ^ 35.5 - 36. Based on the ASCA observations of 
the Small Magellanic Cloud and on the behavior of the in- 
tegrated X-ray luminosity of distant galaxies located at red- 
shifts z ~ 0. 3 - 1.3 observed by Chandra in the Hubble Deep 
Field North. IGrimm et ail -20031) tentatively suggested that the 
HMXB XLF is not dramatically affected by metallicity varia- 
tions. 
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Study of the population of high mass X-ray binaries in the 
LMC is of importance for several reasons: 

1. The Magellanic Clouds ar e known to have a significant 
under-abundance of metals JWesterlundLll997l) . For the in- 
terstellar medium and young stellar population in the LMC 
the metal abundance is approximately ~ 1/3 of the sola r 
value llWesterliincifl997tlCTarnettHl999l:lKorn et all E)02). 
The effects of the metallicity variations on the population of 
X-ray binaries are poorly understood. Study of the LMC, a 
galaxy with a relatively well-known chemical composition 
gives a unique opportunity to investigate such effects obser- 
vationally. 

2. Owing to the proximity of the LMC, the weakest sources 
become reachable with a moderate observing time. Indeed, 
the sensitivity of a typical Chandra or XMM-Newton ob- 
servation, ~ 10" 14 erg/s/cm 2 , corresponds to the luminosity 
of ~ 3 • 10 33 erg/s at the LMC distance. This opens the 
possibility to study the low luminosity end of the HMXB 
XLF, extending by ~ 2-3 orders of magnitud e towards low 
luminosities the luminosity range studied by iGrimm et alJ 
(2003). Among other applica tions, this will allow study of 
the impact of propeller effect (111 arionov & Siinvaevlll975l) 
on the luminosity distribution of HMXBs. 

3. The proximity and moderate inclination angle of the 
Magellanic Clouds (especially that of the LMC) allows 
one to construct Herzsprung-Russel diagrams and accu- 
rately determine ages, star formation history and initial 
mass function of various components of the stellar popu- 
lation. In particular, large aggregates of young and coeval 
stellar populations of ages ranging from £1 — 2 Myr to 
^ 50 - 100 Myr were discovered in the LMC. The most 
extreme examples o f these are the LMC 4 supergiant shell 
jBraun et all ll 997t> and R136 (HD 38 268) stellar cluster 
near the center of the nebula 30 Dor dMassev & Hu nter. 
Il998l) . Combined with HMXB number counts in the opti- 
cally studied regions, these results provide a unique possi- 
bility to directly determine the efficiency of HMXB forma- 
tion as a function of time elapsed since the star formation 
event. 

The proximity of the LMC, on the other hand, creates sev- 
eral difficulties in its studies. Due to its large angular extent 
on the sky, ~ 10° x 10°, multiple mosaiced observations are 
required to achieve a meaningful coverage of the galaxy. The 
low surface density of X-ray binaries results in a large fraction 
of foreground stars and background AGNs among the detected 
sources, thus creating the problem of separating true galactic 
members from interlopers. This was one of the main obsta- 
cles in the earlier studies of X-ray binaries in the LMC with 
moderate angular resolution, insufficient for reliable identifica- 
tion of X-ray sources with optical catalogs. With the advent of 
Chandra and XMM-Newton observatories, the latter problem 
is to a large extent alleviated. 

In the present paper we study the population of X-ray bina- 
ries in the LMC based on archival XMM-Newton data. The dis- 
tance modulus of the LMC is m— M = 18.45+0.10 JWesterlund 
Il997l) . corresponding to the distance of D m 50 kpc. The inter- 
stellar reddening varies across the LMC, with typical values in 
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Fig.l. Far-infrared (IRAS, 100//) map of the 10° x 10° re- 
gion around the Large Magellanic Cloud. Circles show fields 
of view of the XMM-Newton observations used for the analy- 
sis. 

the range E B - V * 0.05 - 0.2 dWesterlundl Il997l) . A value of 
Eb-v = 0.075 corresponding to the direction towards the nom- 
inal center of the LMC is used throughout the paper. 

The paper is structured as follows. XMM-Newton obser- 
vations and data analysis are described in section [3] In sec- 
tion [3]we discuss the nature of detected X-ray sources. In sec- 
tion @]we describe the HMXB search procedure and its results. 
Resulting HMXB luminosity function and the CXB log(N)- 
log(S) are presented in section [5] The impact of the propeller 
effect on the HMXB luminosity function is considered in sec- 
tion [6] In section 0we discuss observed spatial non-uniformity 
of the N hmxb IS FR ratio and its dependence on the age of the 
underlying stellar population. Our results are summarized in 
section|8] 

2. Observations and data analysis. 

We have selected 23 XMM-Newton archival observations with 
the pointing direction towards the LMC and with a sensitivity 
better than ~ few x 10~ 14 erg/s/cm 2 in the 2-10 keV energy 
band. These observations are listed in Tabled Fig[2shows their 
fields of view, overlayed on the far-infrared map (IRAS, 100//) 
of the Large Magellanic Cloud. 

The observations were processed with the standard SAS 
task chain. After filtering out high background intervals we 
extracted images in the 2-8 keV energy band. This energy 
range was chosen to minimize the fraction of supernovae rem- 
nants, cataclysmic variables and foreground stars among de- 
tected sources, generally having softer spectra than high mass 
X-ray binaries. To improve the sensitivity of the survey, images 
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Table 1. List of XMM-Newton observations used for analysis 
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Obs. ID 


Target 


R.A. 


Dec. 


Instrument 


Exposure 


SFR (1) 








J2000 


J2000 




ksec 


M Q /yr 




0113000501 


0519-69.0 


05 19 23 


-69 00 52 


MOS1+MOS2 


25+47 


4.13 


10" 


-3 


0111130201 


2E 0509.5-6734 


05 09 16 


-67 31 41 


MOS1+MOS2 


34+34 


7.77 


io- 


-4 


0062340101 


2E 0453.8-6834 


04 53 41 


-68 31 07 


MOS1+MOS2 


16+16 


1.49 


io- 


-3 


0111130301 


2E 0525.3-6601 


05 25 13 


-65 59 38 


MOS1+MOS2 


9.5+9.5 


1.83 


io- 


-3 


0134520701 


AB Dor 


05 28 35 


-65 28 30 


MOS1+MOS2 


48+48 


2.92 


io- 


-4 


0123510101 


CAL83 


05 43 33 


-68 24 07 


MOS1+MOS2 


10+10 


1.37 


io- 


-3 


0104060101 


LMC Deep field 


05 31 17 


-65 55 58 


MOS1+MOS2 


38+38 


6.35 


io- 


-4 


0109990201 


DEML71 


05 05 49 


-67 54 1 1 


MOS1+MOS2 


22+23 


1.30 


io- 


-3 


0109990101 


LHA 120-N63A 


05 35 51 


-66 00 32 


MOS1+MOS2 


9.6+9.6 


1.38 


io- 


-3 


0094410101 


LMC-2(fl) 


05 42 35 


-69 03 20 


PN 


8.8 


9.74 


io- 


-3 


0094410201 


LMC-2(f2) 


05 42 58 


-69 28 22 


PN 


8.8 


5.51 


io- 


-3 


0094410401 


LMC-2(f4) 


05 46 48 


-69 33 39 


MOS1+MOS2 


5.2+5.2 


4.60 


io- 


-3 


0126500101 


LMC X-3 


05 38 44 


-64 06 1 1 


MOS1+MOS2 


20+20 


7.82 


io- 


-5 


0113000301 


N103B 


05 08 41 


-68 43 53 


MOS2 


23ks 


4.52 


io- 


-3 


0089210701 


N120 


05 18 37 


-69 37 56 


MOS1+MOS2 


36+36 


5.76 


io- 


-3 


0089210901 


N206 


05 31 51 


-70 58 39 


MOS1+MOS2 


24+24 


3.35 


io- 


-3 


0071940101 


N51D 


05 26 05 


-67 27 21 


MOS1+MOS2 


32+32 


3.37 


io- 


-3 


0127720201 


Nova LMC 2000 


05 24 41 


-70 14 00 


MOS1+MOS2 


22+22 


1.34 


io- 


-3 


0113000401 


PSR 0540-69.3 


05 40 02 


-69 21 19 


MOS2 


35 


2.29 


io- 


-2 


0113020201 


PSR J0537-6909 


05 37 57 


-69 08 52 


MOS1 


37 


5.08 


io- 


2 


0008020101 


RXJ0439.8-6809 


04 39 57 


-68 07 27 


MOS1+MOS2 


15+15 


1.83 


io- 


-4 


0089210601 


SNR0450-709 


04 50 26 


-70 48 48 


MOS1+MOS2 


56+57 


6.19 


10" 


-4 


0137160201 


YYMen 


04 58 12 


-75 15 01 


MOS1+MOS2 


87+87 


3.50 


10" 


-4 



1 - star formation rate within the XMM-Newton FOV determined form the FIR flux as described in sectionsl3.3landl3l4l 



from MOS1 and MOS2 detectors were merged. If both MOS 
and PN data were available, we used the data having higher 
sensitivity and smaller number of spurious sources. 

2.1. Source detection. 

The source detection on the extracted images was performed 
with standard SAS tasks eboxdetect and emldetect. The value 
of the threshold likelihood L used in the emldetect task to ac- 
cept or reject detected source was chosen as follows. We sim- 
ulated a number of images with Poisson background counts 
(without sources). Each of the generated images was analyzed 
with the full sequence of source detection procedures using dif- 
ferent values of the emldetect threshold likelihood L and for 
each trial value of L the number of detected "sources" was 
counted. The final value of the threshold likelihood, L = 22, 
was chosen such that the total number of spurious detections 
was $ 3 per 23 images. Note that due to the definition of the 
threshold likelihood in the emldetect task, this value should not 
be interpreted as L — — ln(p), with p equal to the probability of 
detecting a given number of counts due to the statistical fluctu- 
ation of the Poisson noise. 

The obtained images were visually inspected and the 
source lists were manually filtered of spurious sources near 
bright sources and arcs caused by a single reflection. All 
extended sources were removed from the lists. The 2-8 



keV source counts were converted to the 2-10 keV en- 
ergy flux assuming a power law spectrum with photon index 
1.7 and Ntf=6T0 20 crrr 2 . The energy conversion factors are 
ecfMos=2.27-l(L 11 erg/cm 2 and ecfpN=8.0T0~ 12 erg/cm 2 . The 
final merged source list contains w 460 sources. Their flux dis- 
tribution is plotted in Fig|3] The sources with flux > 10 -13 
erg/s/cm 2 are listed in Table|2 

2.2. Boresight correction 

For several XMM observations we performed boresight correc- 
tion using the SAS t ask eposcorr and o ptical sources from the 
USNO-B catalogue jMonet et all 12003b . The correction was 
applied if any of the following conditions was satisfied: 



1 . Optical counterparts for 2 or more well-known sufficiently 
bright X-ray sources were found. 

2. Only one well-known X-ray source with an optical counter- 
part in USNO is present. Offsets calculated by the eposcorr 
task statistically significantly reduce its displacement from 
the optical counterpart. 

3. The offsets calculated by the eposcorr task using optical 
data in different non-overlapping magnitude ranges were 
consistent with each other. 
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Q Q ^ I I I I I ll I I I I I I III I I I I I I 

10 -15 10 -14 1Q -13 10-12 

2—10 keV flux, erg cm -2 s _1 

Fig. 2. The survey area as a function of 2-10 keV flux, calcu- 
lated as described in section l2~3l 

2.3. Correction for incompleteness 

To compute the flux-dependent survey area, we reproduced the 
likelihood calculation procedure used in the emldetect task and 
computed the point source detection sensitivity map for each 
observation. The sensitivity in a given pixel of the image was 
defined as the count rate of the source located in this pixel, 
whose likelihood of detection was equal to the threshold value 
used in the emldetect task. Combining sensitivity maps of all 
observations we calculated the survey area as a function of flux 
(Fig. |2j. From Fig.|2]one can see that the incompleteness effects 
become important at fluxes ^ 10~ I3 erg/s/cm 2 . In the high flux 
limit, the total area of the survey equals A tot ~ 3.77 deg 2 . 

The inverse area of the survey can be used to correct the 
observed differential log(N)-log(S) distribution: 

*») =J-.fa (1) 

dS /corrected MS) \uS ) observed 

where A(S) is the survey area at flux S. The corrected cumula- 
tive log(N)-log(S) distribution can be obtained as follows: 

where S j is the flux of the j-th source. With these corrections, 
the flux distributions are normalized per deg 2 . The flux distri- 
bution, corrected for incompleteness effects, is plotted as the 
thick histogram in Fig|3] 

3. Nature of X-ray sources in the field of LMC 

3.1. Background and foreground sources 

The total number of sources with the flux Fx [2 - lOkev] > 
3 . 34 ■ 1 0~ 1 4 erg/ s/cm 2 is 181. After correction for the survey in- 



completeness (eq|2i this number becomes N b S (> 3.3 ■ 10~ 14 ) « 
214. According to t he CXB log(AT) - log(S) determined by 
iMoretti et alJ J2003I) . the total number of CXB sources ex- 
pected in the field of 3.77 deg 2 is Ncxb ~ 218^^ 7 . The errors 
in this estimate were computed from the uncertainty of the nor- 
malization of the C XB log(A0 - log(5 ) relation, as given by 
IMoretti et alJ J2003h . From the comparison of these numbers it 
is obvious that the majority of the detected sources are back- 
ground AGNs. 

A significantly less important source of contamination 
are X-ray sources associated with foreground stars in the 
Galaxy (no known Galactic X-ray binaries were located in the 
field of view of XMM observations). There were seven well- 
known bright nearby stars with 2-10 keV flux exceeding 10~ 14 
erg/s/cm 2 . All these were excluded from the final source list. 

3.2. Low mass X-ray binaries 

Given the limiting sensitivity of the survey, Fx ~ (1 — 3) • 10~ 14 
erg/s/cm 2 , corresponding to the luminosity Lx ~ 3 ■ 10 33 - 10 34 
erg/s at the LMC distance, the intrinsic LMC sources are dom- 
inated by X-ray binaries. Their total number is proportional to 
the stellar mass (LMXB) and star formation rate (HMXB) of 
the galaxy. 

The total stellar mass of the LMC can be estimated from 
t he integrated optical luminos ity. According to the RC3 catalog 
de Vaucouleurs et all I 199 lh . the reddening-corrected V-band 
magnitude of the LMC equals Vjo * 0.13, corresponding to the 
total V-band luminosity of Ly ~ 1.9 • 10 9 Ln- From the dered- 
dened optical color (B - V)tc, ~ 0-44 Jde Vaucouleurs et all 
1 199 ll) and using results of bell & de Jondd200ll) . the V-band 
mass-to-light ratio is (M/L)y ~ 0.77 in solar units, giv- 
ing the total stellar mass of the LMC M, ~ 1.5 • 10 9 Mq. 
Using more recent determin ation of the (B - V)t ~ 0.54 
(bothun & Thompson! Il988l) . the V-band mass-to-light ratio 
is (M/L)y ~ 1.0, increasing the total stellar mass estimate to 
M, ~ 2.0 ■ 10 9 Mq. For the stellar m ass of * (1.5 - 2.0) ■ 10 9 
M Q , using results of iGilfanovl J2004I) . -8-10 LMXBs with 
luminosity Lx ^ 10 35 erg/s are expected in the entire galaxy. 

3.3. Star formation rate in the LMC 

Below we compare the estimates of the star formation rate in 
the LMC, obtained from different SFR indicators. As is com- 
monly used, the SFR values quoted below refer to the forma- 
tion rate of stars in the 0.1 - 100 M range, assuming Salpeter's 
initial mass function (IMF). 

Based on the total H a luminosity of the LMC and applying 
an extinction correction of A(H a ) » 0.3 Ke nnicutt et alJ (11995) 
estimated the star formation rate of 

S¥R(H a ) = 0.26 M /yr (3) 

As was discussed by iKennicutl ( 1 199 lb . i/^-based SFR indi- 
cators can underestimate the total star formation rate in the 
Magellanic Clouds, due to uncertainty of the escaping fraction 
of the ionizing radiation and an unaccoun ted for contribution 
of the diffuse ionized gas. Kennicutt ( 1991 ) gave an upper limit 
of SFRrS 0.6 M /yr. 
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Table 2. List of detected sources with 2-10 keV flux F x > lCL 13 erg/s/cm 2 
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Source ID 


Fx[2- 


10 keV] 


L x [2-10keV]' 


Rmag 2 




Comments 




erg/ 


s/cm 2 


er 


g/s 








XMMUJ053856.7-640503 


2.67 


■ 10" 


1 


8.0- 


10 36 


16.7 


33 


HMXB LMC X-3 


XMMUJ05401 1.5-691953 


2.45 


■ 10" 


1 


7.4- 


10 36 






Crab-like pulsar PSR B0540-69 


XMMUJ052844.9-652656 


1.30 


■ 10" 


1 






6.4* 


0.0012 


foreground star ABDor 


XMMUJ045818.0-751637 


6.06 


• 10" 


2 






7.5* 


0.0016 


foreground star YYMen 


XMMUJ053747.4-691020 


5.39 


• 10" 


2 


1.6- 


10 36 






Crab-like pulsar PSRJ0537-6910 


XMMUJ053 1 1 3 . 1 -660707 


3.01 


• 10" 


2 


9.0- 


10 35 


14.0 


0.31 


HMXB Be EXO05 3 109-6609 


XMMUJ052402. 1-701 108 


2.08 


• 10" 


2 










AGNRXJ0524.0-7011 


XMMUJ04383 1 .4-68 1 203 


8.88 


• io- 


3 


2.7 • 


10 35 








XMMUJ053844.2-690608 


8.35 


• io- 


3 


2.5- 


10 35 






W-R stars (?) inR136 


XMMUJ05301 1.3-655123 


6.59 


• 10" 


3 


2.0- 


10 35 


14.7 


0.13 


HMXB Be (?) 272s pulsar 3 


XMMUJ050550.2-675018 


5.80 


• 10" 


3 


1.7 • 


10 35 


17.6 


1.5 




XMMUJ053 1 15.4-705350 


5.25 


• 10" 


3 


1.6- 


10 35 


13.8 


0.043 




XMMUJ0532 1 8.6-7 10746 


4.28 


• 10" 


3 


1.3- 


10 35 


17.9 


1.6 




XMMUJ053833.9-691 157 


3.65 


• 10" 


3 


1.1 • 


10 35 


17.8 


1.3 




XMMUJ050026.6-750449 


3.50 


• 10" 


3 


1.1 • 


10 35 


18.2 


1.7 




XMMUJ052847.3-653956 


2.70 


• 10" 


3 


8.1 • 


10 34 








XMMUJ045208.7-705652 


2.33 


• io- 


3 


7.0- 


10 34 








XMMUJ045637.9-751015 


2.06 


• io- 


3 


6.2- 


10 34 


18.3 


1.1 




XMMUJ054 1 34.7-682550 


2.04 


• 10" 


3 


6.1 • 


10 34 


14.0 


0.020 




XMMUJ0523 12.9-701531 


1.94 


• 10" 


3 


5.8- 


10 34 








XMMUJ050526.6-674313 


1.91 


• 10" 


3 






10.8 4 


0.0011 


foreground star GSC 9161.1 103 


XMMUJ053528.5-691614 


1.89 


• io- 


3 


5.7- 


10 34 






SNR SN1987A 


XMMUJ054 1 5 1 .2-682609 


1.75 


• 10" 


3 


5.3- 


10 34 








XMMUJ045244.2-70482 1 


1.71 


• 10" 


3 


5.1 • 


10 34 








XMMUJ050833.2-685427 


1.70 


• 10" 


3 


5.1 • 


10 34 








XMMUJ052403. 1-673606 


1.67 


• 10" 


3 


5.0- 


10 34 


18.3 


0.91 




XMMUJ05304 1.1-660535 


1.66 


• 10" 


3 


5.0- 


10 34 


18.2 


0.83 




XMMUJ053810.0-685658 


1.58 


• 10" 


3 


4.7 • 


10 34 








XMMUJ052947.7-655643 


1.51 


• 10" 


3 


4.5- 


10 34 


14.8 


0.033 


HMXB Be/X transient RXJ0529.8-6556 


XMMUJ0452 1 9. 8-684 1 5 1 


1.49 


• 10" 


3 


4.5- 


10 34 








XMMUJ053734.9-64 10 1 2 


1.46 


• 10" 


3 


4.4- 


10 34 








XMMUJ053909.8-6403 10 


1.33 


• 10" 


3 


4.0- 


10 34 








XMMUJ051832.5-693521 


1.32 


• 10" 


3 


4.0- 


10 34 








XMMUJ05 1747.0-685458 


1.32 


• 10" 


3 


4.0- 


10 34 








XMMU J05 1 003 .9-67 1858 


1.31 


• 10" 


3 


3.9- 


10 34 








XMMUJ054237.8-683205 


1.29 


• 10" 


3 


3.9- 


10 34 








XMMUJ050736.6-684752 


1.23 


• 10" 


3 


3.7- 


10 34 








XMMUJ053257.7-705 113 


1.23 


• 10" 


3 


3.7- 


10 34 


19.0* 


1.3 




XMMUJ052422.4-672020 


1.22 


• 10" 


3 


3.7- 


10 34 








XMMUJ052353.7-672824 


1.21 


• 10" 


3 


3.6- 


10 34 








XMMUJ054212.6-692442 


1.19 


• 10" 


3 


3.6- 


10 34 


19.3* 


1.6 




XMMUJ054543.2-682528 


1.17 


• 10" 


3 


3.5 • 


10 34 








XMMUJ050048.6-750653 


1.15 


• 10" 


3 


3.5 ■ 


10 34 








XMMUJ053632.6-655644 


1.14 


• 10" 


3 


3.4- 


10 34 








XMMUJ045 137.0-710210 


1.14 


• 10" 


3 


3.4- 


10 34 








XMMUJ053841.7-690514 


1.03 


• 10" 


3 


3.1 • 


10 34 






W-R stars (?) in R140 


XMMUJ052558.0-701 107 


1.03 


• 10" 


3 






11.0 


0.00068 


foreground star, K2IV-Vp, RS CVn 


XMMUJ052952.9-705850 


1.02 


• 10" 


3 


3.1 • 


10 34 









1 - calculated assuming distance to LMC 50 kpc. 

2 - apparent magnitude in the R-band of the optical counterpart from GSC2.2.1 catalogue (Morrison & Mc Learl EoOlh . located within 3.6 
arc sec from the X-ray source. 

3 - lHaberl et all 12003 ) 

4 - apparent R-band magnitude from the USNO-B catalogue iMonet et alll2003l) . 
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According to the Catalog of IRAS observations of large op- 
tical galaxies jRice et al.l ll988). the total infra-red luminosity 
of the LMC i s L m « 2.54 • 10 4 2 erg/s. With the IR-based SFR 
calibration o flKennicuttl i 1 998h . this corresponds to the star for- 
mation rate of 

SFR(IR) = 0.23 M /yr (4) 

i.e. it appears to be consistent with the H a -based estimate. 
However, similarly to H a , the IR luminosity luminosity can 
also underestimate the star formation rate in the low-SFR 
galaxies, such as the Magellanic Clouds, due to uncertain and, 
possibly, t he large value of the photon escape fraction (e.g. 
lBdill2003) . An SFR estimator relatively free of this uncertainty 
is the one based on combined IR and UV emission. The total 
integrated U V flux of the LMC wa s measured by the D 2 B- 
Aura satellite ( tVangioni-Flam et allfl980l) . F A = 3.36- 10 9 and 
1.95-10" 9 erg/s/cm 2 /A atA= 1690A and 2200 A respectively. 
Correcting these value for the foreground Galactic extinction, 
(E(B-V)=0 .075,Aifflo = 0.59 a ndA?7nn = 0.72) a nd following 
iBelll J2003I) . we obtain, with the iKennicuttl i 1 998l) calibrations, 

SFR(UV) + SFR(IR) = 0.46 and 0.49 M /yr (5) 

respectively. These values are consistent with the SFR(UV) de- 
rived from the flux es corrected fo r both foreground and internal 
extinction JVangioni-Flam et"al l ll98(t : 

SFR(UV, ext.corr.) = 0.40 - 0.65 M Q /yr (6) 

where the lower and upper li mits correspond to the range o f the 
absor ption corrected fluxe s in lVangioni-Flam et alJ {l]?80). 

lFihpovicetalJ lfl998). from comparison of discrete radio 
(Parkes telescope) and X-ray (ROSAT) sources, estimate a 
number of SNRs in the LMC (36). From the age - radio flux 
density relation they estimate an SNR birth rate of one SNR in 
100 + 20 yrs, corresponding to the star formation rate of: 

SFR(SNR birth rate) * 0.7 + 0.2 M Q /yr (7) 

The SFR values from eq.(|5jQ seem to be more appropriate 
for the physical conditions in the Magellanic Clouds and qual- 
itatively agree with each other. In the following, we assume: 

SFR(LMC) ~ 0.5 + 0.25 M /yr (8) 

The SFR value derived above can be used to predict the ex- 
pected numbe r of high mass X-ray binaries from HMXB-SFR 
calibration of Gri mm et alJ ( 120031) . In applying the HMXB- 
SFR r elation we note that SFR values in the iGrimm et alJ 
(2003) sample are predominantly defined by the FIR and radio- 
based est imators. Th e recent re-calibration of these SFR indi- 
cators by Bell (2003) restored the consistency with other SFR 
indicators and with the observed normalization of IR-radio cor- 
relation, but resulted in the SFR-radio and SFR-IR relation be- 
ing by a factor o f ~ 2.2 and ~ 0.7 different from the com- 
monly used ones JCondoni 1 19921 lKennicutAll998h. W ith this 
new calibration, the SFR values in lGrmimetalT l l2003l) should 
effectively correspond to ~ 1 /3 - 1 /2 of the total formation 
rate of stars in the 0.1 - 100 M mass range. With this in 
m ind, we change the c oefficients in the equations (6) and (7) 
in lGrimm et alJ J2Q 03) to 1.1 and 1.8 respectively and use the 
SFR corresponding to the total mass range, Salpeter IMF, as it 
is commonly assumed. 



3.4. High mass X-ray binaries 

For the star formation rate of SFR=0.5 ± 0.25 Mq yr _1 , we 
predict ~ 60 + 30 HMXBs with luminosities > 10 35 erg s" 1 
in the whole LMC. This confirms that the population of X-ray 
binaries in the LMC is dominated by HMXBs. 

To predict the number of HMXBs in our sample of X- 
ray sources we estimate the star formation rate in the part 
of the LMC covered by X MM pointings from IR AS infrared 
maps provided by SkyView (McGlvnn et all 1 19961) . Calculating 
the far infrared flux according to the formula FIR= 1.26 • 
10 -11 (2.58S6(V + S ^00ll), where FI R is flux in erg/s/cm 2 and 
S a is flux in Jy llHelouetallll985l) . and integrating the IRAS 
maps, we obtained the total SFR= 0.1 1 MQ/yr. This number is 
a factor of ~ 2 smaller than eq.(@}, because it was derived from 
the FIR flux instead of total infrared flux. To make it consistent 
with our determination of the SFR in LMC (eq|SJ, we simply 
multiply the FIR-based values by the correction factor of « 4.5. 
The thus calculated star formation rate of the part of LMC cov- 
ered by XMM observations, excluding r < 4' circle centered 
on R136 (see section l4~5l . is: 

SFR(XMM) * 0.089 + 0.045 M /yr (9) 

Approximately ~ 50% of this value is due to three observations 
whose fields of view included the 30 Doradus giant HII region. 

With the above value of SFR we predict about ~ 1 1 ± 5 
HMXBs with luminosities > 10 35 erg s _1 in the observed part 
of the LMC. The error in this number accounts for the uncer- 
tainty in the SFR estimate and does not include the Poisson 
fluctuations. We note that there already are 5 well-known high 
mass X-ray binaries in our sample: LMC X-3, EXO053109- 
6609, RXJ0529.8-6556, RXJ0532.5-6551 and RXJ0520.5- 
6932 (Tabled- For comparison, the expected number of CXB 
sources with flux F x ^ 3.34 • 10~ 13 erg/s/cm 2 corresponding to 
luminosity L x > 10 35 erg/s, equals N C xb(F x ^ 3.34 ■ 10~ 13 ) * 
6.4. 

There are three "historical" bright high mass X-ray binaries 
in the Large Magellanic Cloud, LMC X-l (L x ~ 1.5 ■ 10 38 
erg/s), LMC X-3 (L x ~ 1.5 • 10 38 erg/s) and LMC X-4 (L x ~ 
0.4 ■ 10 38 erg/s). Their number is consistent with the expected 
value, Nhmxb(Lx > 0.4 ■ 10 38 ) ~ 1.5 + 0.6. Similarly, there is 
one bright LMXB, LMC X-2 with the average luminosity of 
L x ~ 1.5 ■ 10 38 erg/s, which is also consistent with the expected 
number N LM xb(L x > 1.5 ■ 10 38 ) ~ 0.23. 

3.5. Other intrinsic LMC sources 

An important type of X-ray sources associated with star for- 
mation are Wolf-Rayet stars. Due to strong stellar wind with 
M W ind ~ 10 5 — 10 4 MQyr -1 , a W-R star forming a binary sys- 
tem with another W-R or OB star can become a rather lumi- 
nous X-ray source. X-ray emission in such a binary system 
originates from a shock formed b y colliding stellar winds of its 
members (Cherepashchuk, 1976). Typical luminosities usually 
do not exceed L x ^ fewT0 34 erg s _1 . 

There are about 130 known W-R stars in the whole LMC. 
We have cross correlated our data with the fourth cat alogue 
of W-R stars in the LMC by iBrevsacher etaD Jl999l) . Four 
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Fig. 3. Cumulative log(N)-log(S) distribution of detected 
point-like sources, excluding several known foreground stars 
and rotation powered pulsars (section [3}. The thin and thick 
histograms show, respectively, the observed distribution and 
the distribution corrected for the incompleteness effects, as de- 
scribed in section 12.31 The thick gr ey line shows the lo g(N)- 
log(S) of CXB sources according to lMoretti et alJd2Q03l) . 

sources from our sample have been identified as W-R stars. 
Two of them are faint objects with fluxes ~fewT0~ 14 erg 
cm' 2 s _1 identified in optics as WR+OB systems and there- 
fore were excluded from the final source list. The other two 
sources, XMMUJ053844.2-690608 and XMMUJ053833.9- 
691 157, located in 30 Dor, will be discussed in section 1431 

Two well-known Crab-like pulsars located in the field of 
view of XMM-Newton observations of LMC - PSR B0540- 
6910 and PSR J0537-6910- were also excluded from the final 
source list. 

3.6. The Log(N)-Log(S) distribution 

As demonstrated above, the population of compact X-ray 
sources in the LMC field is dominated by two types of sources 
- background AGNs and high mass X-ray binaries in the LMC. 
Their log(N)-log(S) distributions in the flux range of interest 
can be described by a power law with the differential slopes 
* 2.5 (CXB) and * 1.6 (HMXBs). Due to a significant differ- 
ence in the slopes, their relative contributions depend strongly 
on the flux. At large fluxes, Fx ^ (2 - 3) ■ 10~ 13 erg/s/cm 2 
(Lx <: 10 35 erg/s) the X-ray binaries in LMC prevail. On the 
contrary, in the low flux limit, e.g. near the sensitivity limit of 
our survey, Fx ~ 10" 14 erg/s/cm 2 , the vast majority of the X- 
ray sources are background AGN. 

This is illustrated by Fig. showing the observed and 
corrected for incompleteness log(N)-log(S) distribution of all 



sources from the final source list. The corrected log(N)-log(S) 
distribution agrees at low fluxes with that of CXB sources. At 
high fluxes, there is an apparent excess of sources above the 
numbers predicted by the CXB sources log(N)-log(S), due to 
the contribution of HMXBs. 

4. Identification of HMXB candidates 

To filter out contaminating background and foreground 
sources, we use the fact that optical emission from HMXBs is 
dominated by the optical companion, whose properties, such as 
absolute magnitudes and intrinsic colors, are sufficiently well 
known. The X-ray-to-optical flux ratios of HMXBs also occupy 
a rather well-defined range. In addition, we take into account 
the fact that intrinsic LM C objects h ave small proper motions, 
$ 1 - 2 mas/yr ( WesterlundJ Il997ft . which helps to reject a 
number of foreground stars with high proper motion. 

4.1. Optical properties of HMXBs counterparts 

High mass X-ray binaries are powered by accretion of mass 
lost from the massive early-type optical companion. The mech- 
anism of accretion could be connected to either (i) strong stellar 
wind from an OB supergiant (or bright giant) or (ii) an equato- 
rial c i rcumstellar disk around Oe or Be type star (e.g.|Corbet, 
1986t Ivan Paradiis & McClintockl Il994l) . Taking into account 
positions of possible optical counterparts on the Hertzsprung- 
Russel diagra m, the distance mo dulus of the LMC, (m - M)q = 
18.45 ± 0.1 dWesterlundl Il997l) . and foreground and intrinsic 
redde ning towards the LMC, E(B-V)~ 0.1 - 0.2 dWesterlundl 
119971) . we can estimate visual magnitudes of HMXB opti- 
cal companions. OB supergiants and bright giants (luminos- 
ity classes I— II) have absolute magnitudes My ~ — 7 -f -4, 
corresponding to an apparent magnitude in the range my ~ 
11.5 4- 14.5. The position of Oe and Be stars in the H-R di- 
agram is close to the main sequence. The majority of such 
systems have optical companions with a spectral class earlier 
than B3. Indeed, in the catalog of high mass X-ray binaries of 
iLiu et alJ J2000l) . in only two systems was the optical counter- 
part classified as B3Ve and in ~ 4 - 5 - later than B3, out of 
~ 90 HMXBs with known optical counterparts. Therefore the 
absolute magnitudes of majority of Be/X binaries are brighter 
than M v <, -1.6(B3Ve) H- -2.4(B2Ve). This corresponds to a 
visual magnitude my S> 16.0 4- 16.9 and about w 0.1 magnitude 
brighter in the R-band. Accounting for the interstellar extinc- 
tion, Ay ~ 0.3 - 0.6, we conclude that the majority of HMXBs 
have apparent magnitudes brighter than my $ 16.5 - 17.5, 
niR ^ 16.2 - 17.0. Such magnitude filtering rejects the majority 
of AGN typically having fainter optical counterparts. 

As potential optical counterparts of HMXBs belong to 
spectral classes earlier than ~ B3, their intrinsic optical and 
near-infrared colors are constrained by R - I S5 -0.19, B - V S5 
-0.20, V-R ^ -0.08, J-K ^ -0.16, 1-K ^ -0.34. Taking ac- 
count of the interstellar reddening, the apparent intrinsic colors 
will not exceed ^ 0.1 - 0.2. 

The interstellar extinction is known to vary across the 
LMC. In the above estimates, t he up per range of the usu- 
ally quoted values dWesterlundl 1 1997b was used. A several 
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times higher extinction is observed in the 30 Doradus region, 
E(B — V) = 0.3 - 0.6, with the mean value in the central re- 
gion E(B - V) = 0.44 + 0.22, and E(B - V ) = 0.65 towards 
the center of the nebula (Westerlund, 1997). Another compli- 
cation is related to the fact that two (USNO-B and GSC) out of 
the three catalogs used for initial search of optical counterparts 
of X-ray sources (section I4~2l have "holes" with a diameter of 
~ 10' - 15' in the 30 Doradus region. This region is considered 
separately in section ^21 

The above discussion is based on our knowledge of HMXB 
optical counterparts in the Milky Way. The depleted metallic - 
ity in the LMC will, of course, affect the optical properties 
of the HMXB counterparts. From the stellar evolution studies 
it is known that, for LMC metallicity, the effective tempera- 
ture of the early- type stars on the ma in sequence increases by 
~ 0.01-0.05 dex (Schr aer et all ll993). resulting in the intrinsic 
colors being by ~ 0.1 - 0.2 magnitude bluer than the Galactic 
ones. As such, these changes do not affect the efficiency of our 
selection criteria. The difference in metallicity might have a 
more significant effect on formation of HMXB systems. These 
effects have not been studied yet and are one of the subjects 
of this paper. On the other han d, from t he optic al spectros copy 
of 14 known HMXBs in LMC. lNegueruella & CoeT ( l2002l) con- 
cluded that their overall optical properties are not very different 
from the observed Galactic population. This suggests that the 
selection criteria based on the optical properties of the Galactic 
HMXBs would be able to identify most of the HMXBs in the 
LMC, with the exception of a small fraction of peculiar ob- 
jects, which in the case of the Milky Way constitute less than 
~ 5 - 10% of the total population of HMXBs. 

4.2. Catalogs and selection criteria 

We have used the following optical and near-infrared catalogs: 

1. USNO-B. version 1.0 llMonet et alll2003l) 

2. Guide Star Catal og, version 2.2.1 (GSC2.2.1) 

jMorrison & McLeanll200ll) 

3. The CCD survey of the Magellanic Clouds jMassevil2002l) 

4. 2-micron All Sky Survey (2MASS) JCutri et alll2003l) 

5. The point source catalogue towards the Magellanic Clouds 
extracted from the data of the Deep Near-Infrar ed Survey 
of the Southern Sky (DENIS) JCioni et allboOOh 

As a first step, we cross-correlated the XMM-Newton X- 
ray source list with the USNO-B, GSC and lMassevl (2002) op- 
tical catalogs, using a search radius of 3.6 arcsec and the fol- 
lowing selection criteria: R < 17.5,/?-/ < 1.5, B-R < 1.5 and 
V - R < 1.5. If an optical object was present in the USNO-B 
and GSC catalogs, the preference was given to GSC, as it has a 
higher photometric accuracy. The color limits are significantly 
higher than the possible colors of HMXB optical companions 
and were chosen to account for the limited photometric accu- 
racy of optical catalogues. Relaxed color limits also allow for 
variations of the intrinsic LMC reddening up to E(B -V) ~ 1 .5, 
corresponding to N H ~ 8 ■ 10 21 cirT 2 . We have found opti- 
cal counterparts for 78 X-ray sources. Taking into account the 
surface density of the optical objects in the vicinity of X-ray 



sources and the value of the search radius, we estimate the num- 
ber of random matches to be 30, i.e. about half of all optical 
matches. 

In the second stage, we cross-correlated all X-ray sources 
having optical counterparts with the near-infrared catalogs 
2MASS and DENIS, using the same search radius as before. 
The following filtering procedure was then applied. 

1 . For X-ray sources with near-infrared counterparts, all those 
with the infrared colors J-K or I-K >0.7 were excluded 
from the following consideration. 

2. All X-ray sources whose optical counterparts had detected 
proper motion were rejected (all such objects had proper 
motion > 10 mas/yr, significantly exceeding that of the 
LMC). 

3. If the optical counterpart was located at a distance of more 
than two positional errors of the X-ray source plus 1 .5", the 
X-ray sources was rejected. The 1.5" was added to allow 
for the astrometric accuracy of the XMM-Newton bore- 
sight. 

4. All X-ray sources with low X-ray-to-optical flux ratio 
Fx/F opt < 10~ 3 were rejected. The optical flux was calcu- 
lated from the R-band magnitude using F opt = 3.83 ■ 10 -6 ■ 
10 -»W2.5 erg/s/cm 2 . 

Such low Fx/F opt < 10~ 3 ratios are typical for foreground 
stars but not for X-ray binaries. All confirmed HMXBs in 
our sample have Fx/F opt > 10~ 2 , except RXJ0532. 5-6551 
and RXJ0520.5-6932, which have F x /F opt * 3.6 • 10~ 3 and 
« 4.6 • 10~ 3 respectively. 

5. If an optical counterpart was found also in the [jylassev 
J2002li catalogue, having high photometric accuracy, we 
used its optical colors for spectral classification, rejecting 
X-ray sources with counterparts of spectral type later than 
AO. 

4.3. Search radius 

To choose the search radius for cross-correlation with the op- 
tical and infrared catalogs we analyzed the dependence of 
the number of matches between the X-ray source list and 
the USNO-B catalog as a function of the search radius. No 
magnitude or color filtering was performed for this analysis. 
Furthermore, for a noticeable fraction of optical objects, the 
USNO-B catalog contains multiple entries. No attempt to filter 
out such multiple entries was made (unlike for the final analy- 
sis, described in sections l4~2l and I4~4l . Therefore the absolute 
numbers of the total number of matches and of chance coin- 
cidences should not be directly compared with the numbers in 
section l4~2l 

At large values of the search radius, ro ^ 10" - 20", 
the observed number of matches asymptotically approaches 
the A^ ma tch K '"y law, expected for the number of chance co- 
incidences (Fig|4] left panel). As the localization accuracy of 
XMM-Newton is sufficiently high, generally better than ~few 
arcsec, at small values of the search radius, the total number 
of matches is dominated by true optical counterparts of X- 
ray sources. As is obvious from Fig@] the chosen value of the 
search radius, ro = 3.6 arcsec, allows us on one hand to detect 
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Fig. 4. Left: Ratio of the actual number of matches between the X-ray source list and USNO-B catalog to the the expected number 
of chance coincidences, as a function of the search radius. Right: Difference between the number of matches and the expected 
number of chance coincidences. The vertical dashed line in both panels shows the search radius of 3.6 arcsec. The shaded area 
indicates the Poisson uncertainty, computed from the square root of the actual number of matches. 



a significant fraction, ^ 80 - 90%, of true optical counterparts 
(right panel). On the other hand, it results in a reasonable frac- 
tion of chance coincidences, % 50 - 60% (Fig@] left panel). 

4.4. Identification results 

After the filtering, 28 X-ray sources were left in the list of po- 
tential HMXB candidates (Table QJ. Of these, 9 sources have 
optical properties consistent with HMXB counterparts. These 
sources are considered likely HMXB candidates and are listed 
in the first half of the table. All known HMXBs located in 
the filed of view of the XMM-Newton observations are among 
these sources. Listed in the second half are the sources whose 
nature cannot be reliably established based on the available op- 
tical data. 

Comments on the individual sources: 

#9: This source has multiple counterparts in USNO and 
GSC catalogues. One of USNO sources has high proper mo- 
tion and high F v /F 0/ „ ~ 10~ 2 , therefore, it is likely a chance 
coincidence. Another three GSC and one USNO sources from 
significantly different epochs are located at nearly the same po- 
sition and show no evidence of proper motion. Their color in- 
dexes have large uncertainties, but all agree with the HMXB 
nature of the source. It is included in the list of likely HMXB 
candidates, although it is less reliable than other sources in this 
group. 

#12: Th i s sou rce has been previously classified by 
Sasa ki etal] d2000l) as an HMXB candidate because of the 
nearby B2 supergiant. However, its improved X-ray position 
obtained by XMM deviates by « 20" from the proposed B2 
star. This significantly exceeds the positional error of cr r w 2.7" 
making the association with the B2 supergiant improbable. The 



optical matches found within 3.6" do not allow us to reach def- 
inite conclusion about the nature of this source. 

#13: The opti cal counterpart of this source was found in 
the lMassevI J2002 ) catalogue. The color indexes B-V=0.25 and 
V — R = 0.17 and proper motion of the nearby USNO source 
do not allow us to draw conclusions about its nature. 

4.5. 30 Doradus region 

Owing to high stellar density, the central part of this luminous 
HII region is not well represented in the all-sky optical cata- 
logs, such as USNO and GSC. 

There are 5 X-ray sources in t he final XMM list within the 
nominal size of the nebula, r < 7' ( Weste riundlll997() . Of these, 
3 are located outside r « 4.5 arcmin and a search for their op- 
tical counterparts does not present a problem. The remaining 
2 sources, XMMUJ053844.2-690608 and XMMUJ053841.7- 
690514, are located at r < V from the center of 30 Dor 
and positionally coincide with the R136 and R140 stellar 
clusters respectively. Both sources show evidence for non- 
zero angular extent in XMM data. They were detected earlier 
by ROSAT a nd were classified as high mass X-ray binaries 
(|Wan BHl995h . Based on Chandra data, IPortegies Zwart et alJ 
d2002l) resolved XMMUJ053844.2-690608 into one bright 
source and a number of weaker sources, all positionally co- 
incident with bright early type (03f* or WN) stars in the 
R136 cluster. XMMUJ05 3 84 1.7-6905 14 is located in the 
R140 stellar cluster, known to contain at least 2 WN stars. 
Based on the positional coincidence with Wolf-Rayet stars 
and, mainly, on the young age of the R136 cluster, $1-2 
Myr dMassev & Hunte ri Il998l). insu fficient to form compact 
objects, IPortegies Zwart et all J2002I) suggested that all these 
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Table 3. List of HMXB candidates 



# 


R.A. 


Dec. 




(1) 


L x 


(1) 




V-R 




Comments 








erg/s/cm 2 


erg 


s- 1 










Likely HMXB candidates^ 


1 


05 38 56.7 


-64 05 03 


2.67 


lo-' 1 


8.0- 


10 36 


16.7 


-0.2 5 


33 


LMC X-3 


2 


05 31 13.1 


-66 07 07 


3.01 


io- 12 


9.0- 


IO 35 


14.0 


0.23 


0.31 


Be/XEXO053 109-6609 


3 


05 30 11.4 


-65 51 23.4 


6.59 


io- 13 


2.0- 


IO 35 


14.7 


0.18 


0.13 


Be/X ? 272s pulsar (2) 


4 


05 31 15.4 


-70 53 50 


5.25 


io- 13 


1.6- 


IO 35 


13.6 


0.1 


0.038 




5 


05 41 34.7 


-68 25 50 


2.04 


io- 13 


6.1 • 


10 34 


14.0 


-0.03 


0.02 




6 


05 29 47.7 


-65 56 43 


1.51 


io- 13 


4.5 • 


10 34 


14.6 


-0.02 


0.027 


Be/X transient RXJ0529.8-6556 


7 


05 32 32.5 


-65 51 41 


5.93 


io- 14 


1.8- 


10 34 


13.4 


-0.36 


0.0036 


OB RXJ0532.5-6551 


8 


05 20 29.4 


-69 31 56 


4.05 


io- 14 


1.2- 


10 34 


14.1 


0.03 


0.0046 


RXJ0520.5-6932 


9 


05 31 18.2 


-66 07 30 


3.42 


io- 14 


1.0- 


10 34 


16.0 


-1.3 6 


0.023 
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Wolf-Rayet star in R140 



(1) - 2-10 keV band; (2) - lHaberl et al., 2003); (3) - optical counterpart is consistent with being an HMXB; (4) - the optical counterpart 
exists, but the information is insufficient for its reliable classification; (5) - B-V color from Liu et al. 1 2000); (6) - B-V color; 
Comments: USNO-B - counterpart in USNO-B only; GSC - counterpart in GSC only; M - multiple non-overlapping optical and/or infrared 
counterparts; C - ambiguous color indexes; N - non-HMXB nature is likely; 



sources are colliding wind Wolf-Rayet binaries, rather than 
HMXBs. However, for the two brightest sources the X-ray-to- 
bolometric flux ratios, Fx/F\, a i ~ 10~ 5 , exceed by ~ 1-2 
orders of magnitude the typical values for such objects in the 
Gala xy and for other sourc es detected by Chandra in R136 
JPorteeies Zwart et allEooll) . 



Because of this uncertainty we exclude from further consid- 
eration the r < 4' region, centered on the R136 stellar cluster. 
We note that this does not affect our main results and conclu- 
sions. 



4.6. Completeness 

The completeness of the list of HMXB candidates is defined by 
the following factors: 

1. The completeness of the optical catalogues. 

The initial search for optical counterparts is based on 
the GSC2.2 and USNO-B catalogs. The GSC2.2 is a 
magnitude-selected (V < 19.5) subset of the GSC-II cata- 
log (http://www-gsss.stsci.edu/gsc/gsc2/GSC2home.htmi. 
The l atter is complete to J - 21 at high galactic latitudes 
jMorrison & McLeanL EoOlh . The USNO-B 10 catalog is 
believed to be complete down to V = 21 dMonet et all 



P.Shtykovskiy and M.Gilfanov: High Mass X-ray Binaries in the LMC 



11 



HMXBs in LMC 



X 

J 

A 




10 



10 35 1 36 

Lx, erg/s 



10 3 



Fig. 5. The incompleteness-corrected XLF of HMXB candi- 
dates in LMC. The upper histogram shows all sources from 
Table [3] except two in the 30 Dor region; the lower one 
shows likely HMXB candidates (upper part of Table[3}. These 
two histograms provide upper and lower limits for the true 
HMXB XLF. The upper grey line and shaded area show the 
luminosity distribution predicted from the "universal" XLF of 
iGrimm et alJ J2003h extrapolated towards low luminosities and 
its uncertainty. The lower solid and dashed lines show the same 
XLF modified by the "propeller effect" assuming the black hole 
fraction of /bh = and /bh = 0.3 respectively (section^. 



120031) . Completeness of both catalogs is known to break 
down in the crowded regions. One example of such a region 
is the central part of 30 Dor which was excluded from the 
analysis (section l4~5l . As no sensitivity maps for the optical 
catalogs exist, a quantitative estimate of the completeness 
of the initial counterpart search is impossible. However, the 
quoted completeness limits of both catalogs are 2 - 3.5 mag 
better than chosen threshold of 17.5 mag for the optical 
counterpart search. This suggests that the completeness of 
the optical catalogs is unlikely to be the primary limiting 
factor. 

2. The efficiency of the initial search due to statistical and sys- 
tematic uncertainties in the positions of X-ray and optical 
sources. This is probably one of the major limiting factors. 
With the chosen value of the search radius, 3.6", we detect 
~ 80 - 90% of true matches (sectionEOl FiggJ. 

3. The filtering procedure applied to the optical matches found 
in the initial search. This procedure is based on expected 
optical and near-infrared properties of HMXBs and, as dis- 
cussed in section PTT1 it would detect ~ 95% of HMXBs 
listed in the catalog of tLiu et a l. (2000). 

Thus, we estimate the overall completeness of our HMXB 
sample to exceed ^ 70 - 80%. As the latter two factors are 



flux-dependent, they might affect not only the total number of 
detected HMXBs, but also their luminosity distribution in the 
faint flux limit. This effect, however, is not of primary con- 
cern as a much larger uncertainty in the luminosity distribution 
at faint fluxes is associated with the sources of unclear nature 
(FigEJ. 

5. HMXBs and CXB sources in the field of the LMC 

5.1. The luminosity function of HMXB candidates in 
the LMC 

The incompleteness-corrected luminosity distribution of 
HMXB candidates is shown in Fig.|5] The upper and lower his- 
tograms correspond to all sources from Tableland to the likely 
HMXB candidates respectively. These two histograms provide 
upper and lower limits for the true X-ray luminosity function 
of HMXBs in the observed part of the LMC. As is clear from 
Fig-H] they coincide at the luminosity L x ^ (2 4- 3) ■ 10 34 erg/s, 
with the uncertainty of the optical identifications becoming sig- 
nificant only at lower luminosities. We note that the upper his- 
togram steepens at low luminosities where its slope is close to 
that of the CXB sources. Such behaviour indicates that a sig- 
nificant fraction of low luminosity sources of uncertain nature 
might be background AGNs. 

In order to constrain the XLF para meters, we fit the data 
using the maximum likelihood method JCrawford et all fl 970) 
with a power law distribution in the range of luminosities 
Lx > 2.5 • 10 34 erg/s, where both curves coincide. We obtain 
best fit value for the differential slope a = 1.28^ 2^ tne nor ~ 
malization corresponds to N(> 10 35 erg/s) » 5 HMXBs. As 
is evident from Fig|5] the slope of the luminosity distribution 
appears to be somewhat flatter and its normalization smaller 
than predicted from extrapolatio n of the " universal" HMXB 
luminosity function of IGrimm et alJ J2003h - a « 1.6 and 
N(> 10 35 erg/s/cm 2 ) * 11 + 5 (section l3~4l 

However, the XLF flattening is not statistically significant 
in the Lx > 2.5 • 10 34 erg/s luminosity range - the Kolmogorov- 
Smirnov probability for the universal HMXB XLF model is 
~ 40%. Only in the entire luminosity range of Fig. [5] is the 
shape of the luminosity distribution of reliable HMXB can- 
didates (the lower histogram in Fig. [5} inconsistent with the 
extrapolation of the universal HMXB XLF, the Kolmogorov- 
Smirnov test giving the probability of ~ 1%. On the other hand, 
the list of reliable HMXB candidates may be incomplete below 
L x ^ (1 -2)- 10 34 erg/s. 

Due to the limited number of sources and ambiguity of 
their optical identifications at low luminosities, it is premature 
to draw a definite conclusion regarding the precise shape of 
the XLF below L x ^ (1 - 2) ■ 10 34 erg/s and its consistency 
with the extrap olation to the low luminosities of the universal 
HMXB XLF of IGrimm etall d2003l) . We note however that the 
flattening of the luminosity distribution leading to the deficit of 
the low luminosity sources should be expected due to the "pro- 
peller effect". As this effect and its impact on the HMXB XLF 
is of interest on its own, we consider it in detail in section[6] 

Another factor affecting the overall normalization of the 
luminosity distribution - the dependence of the number of 
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Fig. 6. Cumulative logN-logS distribution of CXB sources, ob- 
tained after removal of HMXB candidates. The solid line shows 
the distribution from M oretti et al 

HMXBs on the stellar population age - is discussed in section 

El 

5.2. Log(N)-log(S) distribution of CXB sources 

The log(N)-log(S) distribution of CXB sources obtained after 
removal of HMXB candidates is shown in Fig|5] The difference 
between the upper and lower limit on the HMXB population is 
insignificant in this context, due to the large number of CXB 
sources at low fluxes. Below we give values obtained after re- 
moving all HMXB candidates. 

We fit the resulting log(N)-log(S) distribution in the flux 
range F x > 2 ■ 10~ 14 erg cirT 2 s _1 with a power law model 
N(> S) = k(S/S Q )- a , where S Q = 2 ■ 10~ 14 erg/s/cm 2 . Best fit 
values are a = 1.62 ± 0.08 and k = 127 ± 11. The Kolmogorov- 
Smirnov test accepts this model, giving a K-S probability of 
~ 16%. These best fit values agree with the CXB parame- 
ters determined in o ther surveys, in particular with those from 
iMoretti et alJ J2003h : a = 1.57+°;^, k=12r®. The latter dis- 
tribution is shown in Fig|6]by the solid line. 

6. Propeller effect and HMXB XLF 

6.1. The "propeller effect" 

As high mass X-ray binaries are young objects, the neutron star 
magnetic field is sufficiently strong to be dynamically impor- 
tant in the vicinity of the neutron star. Indeed, the majority of 
known HMXBs in the Milky Way Jliu et alll2000l) and Small 
Magellanic Cloud JCorbet et alll20 04) are X-ray pulsars. In the 
presently accepted picture the transition from disk-like accre- 
tion to a magnetospheric flow, co-rotating with the neutron star, 



occurs in a narrow region located at the magnetospheric radius 
R m . The location of the transition region is defined by the bal- 
ance between the NS magnetic field pressure and the pressure 
(ram and thermal) of the accreting matter. There is some uncer- 
tainty in the definition of R„„ due to uncertainty in the physics 
of the disk-magnetosphere interaction, the canonical value be- 
ing (e.g. lLamb et allll973l) : 

R m = I A- ^R^M^L^ cm (10) 

where R(, is the neutron star (NS) radius in units of 10 6 cm, M\ .4 
is its mass divided by \AMq, B\2 is strength of the magnetic 
filed on the NS surface in units of 10 12 Gauss and L35 is X-ray 
luminosity in units of 10 35 erg/s. It was assumed above that the 
mass accretion rate M is related to the X-ray luminosity via 
L X = (GM NS /R NS )M. 

The character of the disk-magnetosphere interaction de- 
pends critically on the fastness parameter, a> = Q*/Qx(i? m ), 
defined as the ratio of the neutron star spin frequency £2, and 
the Keplerian fr equency at the ma gnetospheric radius Q^(7? m ). 
As suggested by Illario nov & Sunvaevl ( 11975ft . at low mass ac- 
cretion rates, the spin frequency of the neutron star can exceed 
the Keplerian frequency at the magnetospheric radius. In this 
case, corresponding to a> > 1, the flow of the matter towards 
the neutron star will be inhibited by the centrifugal force ex- 
erted by the rotating magnetosphere and the matter can be ex- 
pelled from the system due to the "propeller effect". 

If the critical value of the fastness parameter u) prop , at which 
the "propeller effect" occurs, is known, the corresponding value 
of the critical luminosity can be computed. Using ea.( ll0> : 

L x , prop = 3.4 ■ 10 33 a>-™& M-J 3 B 2 12 P^ 3 erg/s (11) 

where P100 is the NS spin period in units of 100 sec. The value 
of Lx, pr0 p defines the lower limit on the possible X-ray lumi- 
nosity of an X-ray binary with given parameters of the neutron 
star. At lower values of M the system will not be observable 
as an X-ray source. The details of the disk-magnetosphere in- 
teraction are not well understood, especially at large values of 
the fastness parameter. It is usually assumed that the "propeller 
effect" occurs at a> > I. However, for the accreting matter to 
be expelled from the system, the linear velocity of the rotat- 
ing magnetosphere should at least exceed the escap e velocity, 
v psr - y/2vK, which corresponds to a> > V2. Spruit & Taam 
Jl993h argued that, due to gas pressure effects, accretion is still 
possible in a narrow range to S 1 . They suggested that the crit- 
ical value of the fastness parameter is a> ^ ylvxlc s , where c s is 
the speed of so und at the magnetospheric r adius. For a standard 
disk solution JShakura & Sunvaevl Il973l) . this value can be 
very high, a> » 1. On the other hand, the standard disk solu- 
tion might be inapplicable i n the v icinity of the magnetospheric 
radius (e.g. S pruit & Taarrl fl993l) . In addition, the disk might 
be rather hot at the inner edge, due to d issipation of the shock 
induc ed by the rotating magnetosphere Jlllarionov & Sunvaevl 
1975). Owing to the strong dependence of the luminosity on 
the fastness parameter (eg. I ll 11 ). this uncertainty in the criti- 
cal value of a> translates into a much greater uncertainty in the 
value of the critical luminosity. With this in mind, we study 
qualitatively the impact of the propeller effect on the luminos- 
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Fig. 7. Distribution of spin periods (left) and magnetic field strength (right. ICoburn et all (2002)) in high mass X-ray binaries. 
The spin period distributions are shown separately for the Milky Way ( iLiu et al.l d200(j) N ) and SMC JCorbet et alJ J2004I) ') sources. 
The SMC distribution is shifted along the x axis by 0.2 dex for clarity. The smooth solid lines in both panels show the model 
distributions used for calculation of the propeller effect. 



ity distribution of X-ray binaries with the neutron star primary, 
assuming that the propeller regime occurs at a> > y/2. 



6.2. Impact on the HMXB luminosity distribution 

The existence of the lower limit on the luminosity of an accret- 
ing neutron star will obviously modify the shape of the lumi- 
nosity function of high mass X-ray binaries, leading to a deficit 
of low luminosity sources. The modified luminosity distribu- 
tion can be calculated as: 

dN = (dN_\ I _R NS 
dL 



dM \GM, 



INS 



f(L) 



dN\ 
dL) 



(12) 



where the luminosity-dependent factor f(L) accounts for 
the "propeller effect". The undisturbed luminosity function 
(dN/dL\) is defined by the the distribution of the binary sys- 
tems over the mass transfer rate M and depends on the distri- 
butions of the binary system parameters and the parameters of 
the optical companions in HMXBs. The function f(L) is given 
by: 



f(L) 



II 



dn dn 



®(L - L pro JP, B)) dPdB 



(13) 



dP dB 

where 0(x) is the Heavy side step function and dn/dP and 
dnjdB are distributions of the HMXBs over the NS spin peri- 
ods and surface magnetic fields, normalized to unity. The min- 
imal luminosity of HMXB L prop {P, B) is given by eq. dl li . with 
a) = V2 and assuming for simplicity that all neutron stars have 
the same mass and radius. 

For the undisturbed distribution (dN/dL) we extrapolated 
towards low luminosities the universal luminosity function of 



HMXBs, derived bv lGrimm etaD J2003h : 



dN 

Tl 



- AL 



■1.6 



(14) 



As it is mainly based on the luminosity distribution of high 
luminosity systems, log(Lx) ~ 36-40.5, it should be relatively 
unaffected by the "propeller effect" and provides a reasonable 
approximation of the M distribution in HMXB systems. 

To estimate the distributions of HMXBs over the NS 
spin period and magnetic field, we used the data on the 
known HMXBs. Strictly speaking, the observed distributions 
are themselves modified by the "propeller effect", as the crit- 
ical luminosity depends strongly on the NS period and mag- 
netic field. For the purpose of this qualitative consideration we 
ignore this effect and use observed distributions. For the spin 
frequency we used the measured period s of known X-ray pul- 
sars in the Milky Way (Liu et al., 2000) and Small Magellanic 
Cloud dCorbet et alll2004l) . The distributions are plotted in the 
left panel of Fig [7] demonstrating that they are similar. We ap- 
proximated these distributions with an empirical function 



dn 
~dP 



1+P/5 + (P/300) 4 



(15) 



where P is the spin period in seconds and the constant «o is 
defined to normalize the distribution to unity. This approxi- 
mation is shown in Fig Q by the solid line. To estimate the 
dnjdB distribution we use the results of the determination of 
the magn etic field strength f rom observations of the cyclotron 
lines bv lCoburn etaD J2002I) . The distribution is shown in the 
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right panel in FigQand was approximated by the log-Gaussian: 



dn 
dB 



no exp 



f(log(B)-log(B )) 



2\ 



2o% 



(16) 



#min < B < B 

n 



with parameters log(Bo) = 12.4 and erg = 0.2. 

The examples of the luminosity distributions, modified by 
the "propeller effect" for different values of the parameters 
are shown in FigFJ] Given the shapes of dn/dP and dnjdB 
distributions, the function f(L) and resulting luminosity func- 
tion weakly depend on the values of Z? max ^ 10 13 Gauss and 
^min 55 0.5 sec. Therefore they were fixed at these values. The 
dependence on P max and B m [ n is significantly stronger. The long 
period pulsars with P <; 10 3 sec are unaffected by the "propeller 
effect" but they can be subject to a significant observational 
bias, as the long periods are more difficult to detect. Based on 
the observed period distribution of the known X-ray pulsars, 
Anax ~ (1 _ 2) • 10 3 sec seems to be a reasonable choice. There 
are no X-ray pulsars with measured B below ~ 10 12 Gauss. On 
the other hand, for B ~ 10 11 - 10 12 Gauss, the cyclotron line 
energy is in the ~ 1 — 10 keV energy range, where it could have 
been easily detected by numerous experiments, operating in the 
standard X-ray band. We assumed in the following B m \„ — 10 12 
Gauss. 

Although the choice of P max and B min significantly affects 
the global shape of the luminosity function, it has a modest ef- 
fect in the luminosity range of interest, \og{Lx) <: 33.5 (Fig|SJ. 
In the simplified model considered above, the shape of the lu- 
minosity function in this range is sensitive only to the radius of 
the neutron star and to the critical value of the fastness param- 
eter u prop (eq.dllj)- 

6.3. Comparison with the observed XLF 

The overall impact of the "propeller effect" on the HMXB lu- 
minosity function depends on the fraction of neutron star bi- 
naries, as HMXBs with a black hole primary are, obviously, 
unaffected by it. Presently, there is some ambiguity in this frac- 
tion. In the Milky Wa y, for 50 out of 8 5 HMXBs, X-ray pul- 
sations were detected (I Liu et all [2000 ). For only a few of the 
remaining 35 sources was the black hole nature of the primary 
confirmed. Therefore, the fraction of HMXBs with black hole 
primaries is 0.04 ^ fbh,MW ^ 0.45. In the case of the Small 
Magellanic Cloud, there are 25 objects listed in the catalog of 
lLiuetal.1 d2000t), of which for 18 X-ray p ulsations were de- 
tected JLiu et all l2000t ICorbet et all 120041) and for no source 
was the black hole nature of the primary confirmed, therefore 

fbh,SMC ~ 0.28. 

With this uncertainty in mind, we compare the model with 
the observed luminosity function of HMXBs in the LMC in 
FigE The luminosity distributions with account taken of the 
"propeller effect" are shown as solid and dashed thick grey 
lines, computed assuming /bh = and /bh = 0.3 respectively. 
The "propeller effect" in both curves was computed assuming 
M m = 1.4 M G , Rns = 15 km, B min = 10 12 Gauss, B max = 10 13 
Gauss, P min = 0.5 sec, P max = 10 3 sec. The f B H - curve illus- 
trates the maximum amplitude of the impact of the "propeller 




Fig. 8. Impact of the propeller effect on the luminosity distri- 
bution of high mass X-ray binaries with a neutron star primary. 
The straight line shows the initial power law distribution with 
differential slope 1.6. The curved lines show the luminosity 
function modified by the propeller effect for different param- 
eters of the distributions of the neutron star spin and magnetic 
field. The neutron star mass and radius are M^s = 1 .4M and 
Rns = 13 km. The dashed line corresponds to our sensitivity 
limit. 

effect" on the luminosity distribution of HMXBs (for the given 
choice of the NS parameters and u) pr op)- 

As expected, the "propeller effect" results in the deficit of 
low luminosity sources and flattening of the XLF. This be- 
haviour is qualitatively similar to the observed XLF (Fig|5}. 
However, due to lack of distinct features of the "propeller ef- 
fect" in the XLF at log(Lx) ^ 33 and large uncertainty in 
the observed HMXB XLF at low luminosities, it is prema- 
ture to draw any definite conclusion regarding the relevance 
of the "propeller regime" of accretion onto a strongly magne- 
tized neutron star. Its impact would be more apparent in the 
\og(Lx) ^ 32 luminosity range. In order to construct this, a 
larger number of sources and better sensitivity limits are re- 
quired. 

7. HMXBs and the age of underlying stellar 
population 

7.1. Spatial distribution of HMXBs in LMC 

Considered globally, the number and luminosity distribution 
of HMXBs in the observed part of the LMC roughly agree 
with expectation based on the universal luminosity function of 
HMXBs, probably modified by the propeller effect at the low 
luminosity end (Fig[5j. However, examination of the numbers 
of detected HMXBs in the individual XMM pointings reveals 
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a significant non-uniformity of the Nhmxb/S FR ratio, as il- 
lustrated by Fig|5] This figure also shows all known HMXBs 
and HMXB candidates in LMC. These do not represent a flux- 
limited sample, therefore, it should be interpreted with caution. 
Nevertheless, their spatial distribution shows the same trends 
as the distribution of the HMXBs detected by XMM, which do 
constitute a flux-limited sample. Fig |9] reveals the remarkable 
absence of correlation between the surface density of HMXBs 
and the star formation rate, as traced by H a emission. Indeed, 
nearly half of the known HMXBs are located in the north-east 
part of the supergiant shell LMC 4 and beyond it. On the other 
hand, there are only a few HMXBs in and near the 30 Dor re- 
gion, which has the highest surface density of star formation, 
as traced by H a and FIR emission, and accounts for w 50% of 
the total SFR in the part of the LMC covered by the survey. 
Factors such as obscuration and the SN kicks do affect the spa- 
tial distribution of HMXBs, but they seem to be insufficient to 
explain the observed distribution. Indeed, the largest interstel- 
lar extinction in the LMC regi on is observed tow ards the center 
of 30 Dor, E(B - V) = 0.65 JWesterlundLll997l) . correspond- 
ing to the hydrogen column density of Nh ~ 3.5 • 10 21 cirT 2 . 
This value cannot affect the spatial distribution of the sources 
in the 2-8 keV band in any significant way. The typical SN 
kick velocities of HMXBs are unlikely to exceed ~ 50 km/s, 
which corresponds to ~ 3.5 arcmin/Myr at the LMC distance. 
Such a velocity is obviously insufficient to significantly modify 
the HMXB spatial distribution on the angular scale of degrees 
over a few million years. A more attractive explanation is of- 
fered by the effects related to the age of the underlying stellar 
population, as discussed below. 

7.2. HMXBs and the stellar age 

HMXBs are often considered as an "instantaneous" tracer of 
star formation. However, the lifetime of the most massive stars 
of M ~ 125 M is ~ 2.6 Myr and can be as long as ~ 15 - 20 
Myr for M ~ 8-10 M Q stars. To this time one should add a de- 
lay required for the binary to reach the HMXB phase, Tps n ~ 1 
Myr. Following an instantaneous creation of stars at t = 0, the 
number of HMXBs Nhmxb(j) is a non-monotonic function of 
time r passed since the star formation event. It equals zero at 
r « 1— 3 Myr, unless stars are formed significantly heavier than 
the conventionally accepted upper mass limit of ~ 125 M Q . At 
later times, the Nhmxb(j) is an increasing function of r, until 
at least r ~ 20 Myr, corresponding to the lifetime of the least 
massive stars, M ~ 8 Mq, capable of leaving behind a compact 
object. The behaviour of Nhmxb{t) at later times is not clear. 
One might expect that at r ^ 20 Myr, the number of HMXBs 
will decrease with time. However, this behaviour will be af- 
fected by the binary systems with a less massive companion, 
entering the HMXB phase at later times (i.e. having a longer 
time delay tpsn)- 

This simple picture might qualitatively explain the ob- 
served non-uniformity of the spatial distribution of HMXBs in 
the LMC. Indeed, the central region of 30 Dor has a very young 
age, m 1 - 2 Myr, a s revealed by its H-R dia gram and stel- 
lar mass distribution (Mass ev"& HunteJ. fl9 98). This is insuf- 
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Fig. 9. The H a map of LMC (greyscale image). The overlayed 
symbol s are: "+"(mage nta) - all known HMXBs from the cat- 
alog of iLiu et alJ (l200Cf> a nd HMXB candidates discovered by 
ROSAT (Kah abkatbOPl : "x" (blue) - likely HMXB candi- 
dates (upper part of Table [3}- Large circles indicate fields of 
view of XMM observations. The oval-shaped region in the 
north-east shows the location of the LMC 4 supergiant shell. 
The region of the highest surface brightness, located to the 
south of Dec= -69°, is the 30 Doradus HII region. 

ficient to form compact objects - neutron st ars or black holes, 
as was noted by Portegi es Zwart et alJ ( 120021) . unless stars with 
masses exceeding ~ 150 M were formed in R136. The LMC4 
supergiant s hell, o n the other hand, has an age of ~ 10 Myr 
(Bra un et al.l ll997). which is sufficient for all stars heavier than 
~ 20 Mq to have become collapsed objects. The OB associa- 
tions in the southern part of LMC 4 are somew hat younger , 
with ages ranging from ~ 2 Myr to ~ 6 - 9 Myr llBraun et all 
1997, and references therein), which might explain the pausity 
of HMXBs there. 

In the b o w- sho ck induced star formation model of 
Ide Boer et al.l dl998h . the stellar population age increases 
clockwise from 1 Myr to the south-east of 30 Dor, to 
~ 10-15 Myr in the north-east of the LMC, near LMC 4, to 
~ 20 - 40 Myr to the north-west. This also qualitatively agrees 
with the observed distribution of HMXBs in the LMC. 

7.3. Number of HMXBs as a function of time after the 
star formation event 

The proximity of LMC allows to construct H-R diagrams and 
sufficiently accurately determine the age and the mass func- 
tion of the stellar population in various parts of LMC. Due 
to the small inclination angle and, consequently, small depth, 
the individual stellar associations can be studied without sig- 
nificant projection effects and contamination by the foreground 
and background populations. As a result of these studies, it has 
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Fig. 10. The time dependence of the specific number of 
HMXBs, rjHMXBij)- The point with error bars and the upper 
limit are data on the LMC 4 supergiant shell and R136 star 
cluster (Table 0}. The solid line shows the simple model de- 
scribed in section 1731 with the normalization deter mined from 
the NhmxbSFR relation from Gri mm et alJ (2003). 

been found that the supergiant shells in the LMC often host co- 
eval stellar populations. One of the most intriguing examples is 
the LMC 4 supergiant shell, with a diameter of ~ 1.4 kpc, in- 
side which no significant age gradients have been found, with 
all th e stars having a pproximately the same age of ~ 9-12 Myr 
jBraun et all fl997). X-ray observations and HMXB number- 
counts in the fields, well studied in the optical band, open a 
unique possibility to directly determine the number of high 
mass X-ray binaries as a function of time elapsed since the star 
formation event. This possibility is explored below using two 
stellar associations in the LMC as an example - the R136 stel- 
lar cluster and the northern part of the LMC 4 supergiant shell. 

We assume that stars are formed instantaneously at time 
t = with the Salpeter IMF with an upper mass cut-off M u 
1 25 Mq. The time-dependent specific number of HMXBs, 
t]hmxb{t), is defined as the number of HMXBs present after 
time t since the star formation event, normalized to the total 
mass of massive stars with M > 8 Mq, formed at t = 0: 

mMxBiT) = m(>Z q b{ !t=o) (17) 

The age of the stellar populatio n in the northern part of LMC 4 
was determined by Bra un et al.l (11997) from the analysis of the 
H-R diagram, t g 10 - 12 Myr, w ith little dispersion between 
individual fields. lBraun et al.|lll997l) also determined the stellar 
mass function and have shown that its slope in the ~ 3 — 15 Mq 
range is close to the Salpeter value of 2.35. The surface mass 
density of the initially-formed massive stars with masses be- 
tween 8 and 125 Mq is w 85 Mq /arcmin 2 . Extrapolating these 



results to the XMM observation of the LMC Deep field, cov- 
ering ss 600 arcmi n 2 and whose FOV did not coincide with 
the fields studied by Br aun et all dl997l) but was almost within 
the boundaries of LMC 4, we can estimate the total mass of 
massive stars within the XMM FOV, formed * 10-12 Myr 
ago, M(> SMq , t = 0) ~ 5 ■ 10 4 Mq. In this observation, 5 
HMXB candidates with luminosity Lx > 10 34 erg/s were de- 
tected. With these numbers we estimate 



t]hmxb(t = 10 - 12Myr) ~ 1 ■ 10~ 4 HMXB/M . 



(18) 



A similar estimate can be made for the R136 regi on. For the 
R136 c luster we complimented the resul ts oflMassev & Hunter! 
(1998) with the flanking fields data of Siria nni et alJ ([2000), 
covering an area of 4.9 arcmin 2 centered on the R136 stellar 
cluster. In this region only two X-ray sources were detected by 
XMM, as discussed in section l4~5l As their nature is uncertain, 
we assume conservatively that the total number of HMXBs in 
this region is < 2. The results are summarized in Tableland 
plotted in FigfTol 

In a simple ad hoc model of the formation of HMXB s, the 
depend ence t]hmxb(t) can be estim ated as follows.lBraun et alJ 
d!997l) . fitting the stellar tracks of Schrae r et alJ 11993), have 
found that the lifetime of a star of mass M is 



U ~ 86 AT 72 Myr 



(19) 



where the mass of the star M is expressed in solar units. For a 
Salpeter mass function, the rate of type II supernovae depends 
on the time t elapsed since the star formation event as follows: 

dN SN 



oc T 88 , 2.6 < t < 19 Myr. 



dt 



(20) 



In the latter inequality t = 2.6 and 19 Myr are the lifetimes 
of the most massive stars and of the stars with M — 8 Mq 
respectively. Taking into account that the HMXB lifetime th is 
short, we obtain for the number of HMXBs active at time r: 



m i ^ dN ™ , x 
Nhmxb(t) — : — (t - T PSN ) ■ t h 

dt 



(21) 



where tpsn is the post-supernova time required for the binary 
to reach the HMXB phase. Here we have assumed that Tpsn 
and th are the same for all binaries. With ea. (l20> we obtain: 



i]hmxb(t) k (t - tpsn) 



I). NX 



(22) 



This equation is valid for 2.6 < t - t psn < 19 Myr. The be- 
haviour of t]hmxb(j) outside this interval, especially at larger t, 
is unclear, as discussed in section l7~2l The shape of t]hmxb(j) 
calculated using eq.|22]with tpsn - 1 Myr is shown in FigllOl 
The normalization in equation i22\ depends on a number of 
parameters of the binary evolution, such as the total fraction of 
the binaries, the fraction, survived the supernova explosion etc, 
whose values are unknown. On the other hand, normalisation 
can be determin ed from t he calibration of the A^mzb-SFR rela- 
tion obtained by Grimm et al . (2003). From its derivation, this 
calibration corresponds to the case of a steady star formation 
on a time scale longer than HMXB formation and life times. 
Therefore: 



Nhmxb 
SFR 



I 



t]hmxb(t) dr. 



(23) 
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Table 4. Constraints on t] HM xb(t) 



8. Summary 



Name 


Age, Myr 


Mi l> 


N {2> 

HMXB 


^HMXB 


R136 


a 1 - 2 


2 ■ 10 4 


< 2 


< 1 • lO 4 


LMC 4 


a 10- 12 


5-10 4 


5 


(1 ±0.45)- 10~ 4 



(1)- Total mass of massive stars, M > 8 Mq; (2) - number of HMXB s 
detected by XMM; (3) - number of HMXBs per unit s t ellar m ass; 
References for age and IMF: R136 - |Massey& Hunter] J 19981) and 
ISirianni et all OoOO), LMC 4 - iBraun et alHl997h 
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Fig. 11. The log(N)-log(S) distribution of X-ray sources in the 
LMC field. The difFerent histograms show the distribution of 
all sources, high mass X-ray binaries and CXB sources. The 
shaded area in the distribution of HMXBs corresponds to the 
range of uncertainties discussed in section lBTl 



The normalization of the curve shown in Fig^H was ob- 
tained from integration of ea. (l23> with t]hmxb(t) defined by 
eq. J22> and the integration limits 2.6 and 19 Myr. The ratio 
Nhmxb/SFR = N HM xb{> 1 34 erg/s)/SF R « 1.9 • 10 3 was 
calculated using eq.(7) from Gri mm et alJ ( 120031) for HMXBs 
with Lx > 10 34 erg/s. As can be seen from Fig EH within the 
accuracy of our (very crude) approximation, there is a good 
agreement with the data. 

The above model is of course very simplified. In its more 
realistic version the effects of spread in Tpsu, thmxb and 
in ages of the stellar population should be taken into ac- 
count. However, to a first approximation such effects will only 
smooth the edges, leaving the overall behaviour of t]hmxb(t) 
unchanged. Significantly more important are the effects of evo- 
lution of the secondary and contribution of the intermediate 
mass systems. 



Based on the archival data of XMM-Newton observations, we 
studied the population of compact sources in the field of the 
LMC. The total area of the survey is w 3.8 sq.degr. with a lim- 
iting sensitivity of ~ 10~ 14 erg/s/cm 2 (Fig®, corresponding to 
the luminosity of ~ 3 • 10 33 erg/s at the LMC distance. 

1. Out of the 460 compact sources detected in the 2-8 keV 
energy band the vast majority, ^ 94%, are CXB sources, 
observed through the LMC (section lXTl Figs Oll it . 

2. Based on the stellar mass and the star formation rate of the 
LMC we demonstrate that the majority, ^ 75 - 90%, of 
the intrinsic LMC sources detected in the 2-8 keV band are 
high mass X-ray binaries (section®. 

3. The proximity of the LMC and the adequate angular res- 
olution of XMM-Newton make it possible to reliably fil- 
ter out the sources whose properties are inconsistent with 
being an HMXB. Based on the optical and infrared mag- 
nitudes and colors of the optical counterparts of the X- 
ray sources we identify 9 likely HMXB candidates (6 of 
which were previously known HMXBs or HMXB candi- 
dates) and 19 sources of uncertain nature (section® Table 
|3|l. The remaining ~ 440 sources, with a few exceptions, are 
background objects, constituting the resolved part of CXB. 
Their flux distribution is consistent with other determina- 
tions of the CXB log(N)-log(S) (section^ Fig®. 

4. With these results we constrain lower and upper bounds of 
the luminosity distribution of HMXBs in the observed part 
of the LMC. We compare these with the extrapolation to- 
wards low luminosities of th e universal lu minosity func- 
tion of HMXBs, derived by iGrimm et all J2003I) . At the 
high luminosity end, the observed distribution is consistent 
with the extrapolation of the universal XLF. The number of 
bright, L x > 10 38 erg/s, "historical" HMXBs in the LMC 
also agrees with the value predicted from its global star for- 
mation rate. At lower luminosities, a deficit of sources is 
observed - the luminosity distribution seems to be some- 
what flatter than the universal XLF (section l5~Tl Fig®. 

5. We consider the impact of the "propeller effect" on the 
luminosity distribution of HMXBs (section® Fig® and 
demonstrate that it can explain qualitatively the observed 
flattening of the XLF at low luminosities (Fig®. We note, 
however, that due to the relatively small number of HMXB 
candidates in the observed part of the LMC and the lack of 
distinct signatures of the "propeller effect" in the \g(Lx) ^ 
33 luminosity domain (Fig®, it is premature to draw a 
definite conclusion about the relevance of the "propeller 
regime" to the accretion onto a strongly magnetized neu- 
tron star. 

6. We found significant field-to-field variations in the number 
of HMXBs, which appear to be uncorrelated with the star 
formation rates inferred by FIR and H a emission (Fig®. 
We suggest that these variations are related to different 
ages of the underlying stellar population in different XMM- 
Newton fields. Using the existence of large coeval stellar 
aggregates in the LMC we constrain the number of HMXBs 
per unit stellar mass as a function of time elapsed since the 
star formation event (section® FigllOi. Based on a sim- 
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pie ad hoc model, we obtain the theoretical dependence 
J]hmxb( t ) an d show that its normalization, as constrained 
by the LMC data, is consistent with t he calibration of the 
NhmxbSFR relation, derived bv lGrimm et al.l (12003). 
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